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1. INTRODUCTION 


IN a former paper (Proc. Ind. Acad. Sci., 1952, 35 A, referred to hereafter 
as p) we considered the following problem: Given that the set of variables 
...-, Xn) Satisfies the conditions: 


(a) 0< eee Xn 


where k,, ....kn_, are fixed but unspecified constants. 


(c) Maximum and minimum values of the successive x’s 
occur alternately, i.e., if x, has its maximum value, x, has its minimum 
value, x; its maximum value and so on. 

To find, under these conditions, » numbers (a, a2...., a,) such that 


It was shown in p that the a’s are the roots of certain Tschebyscheff poly- 
nomials and that they possess certain symmetry properties. 


The object of the present note is to offer some supplementary remarks 
indicating an alternative formulation of the problem as a problem of condi- 
tional maxima and minima. This formulation shows that the original problem 
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is but one of a class of similar problems, all of which are governed by the 
same extremising condition and that the “two extreme solutions” copn- 
sidered in the previous paper (p, p. 212) are indeed such. We give solutions 
for some of these other problems and also present some additional considera- 
tions regarding the original problem. 


2. NEW FORMULATION OF THE PROBLEM 
Given that 


ag< a4< Ag... <= an 


To make ¢ = X,X2....X,= ky an extremum. 


We assume that no two of the a’s in (1) can be equal, for if two of the 
a’s are equal, say a, = a,_,, then x, would be reduced to a constant and 
the problem would be reduced to one in m — 1 variables. 


By the Lagrange method of undetermined multipliers it is readily seen 
that the stationary values of ¢ occur when 


IT (xp — Xs) =0 (3) 


i.e., extreme values of ¢ occur when any two of the x’s are equal and, on 
account of (1), this can happen only when 


Xr = Xr, p= 2, (4) 


Thus ¢ has n — | stationary values. It is now evident that our problem 


in p is that particular case corresponding to the maximum number of equal 
pairs of roots, viz., 


(i) Xy = Xo => O15 X3 = X4 = a3, etc. 
or (ii) Xy = ao, Xo = X3 = Ao, X4 = X5 = 4, etc. 


so that ¢ has only two stationary values. From (3) we see that ¢ is stationary 
when the polynomial 


(x) = x” — kyx™1 + — kgx™3 4 hin gx F kn 
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has a pair of equal roots, the condition for which is the vanishing of the elimi- 
nant of 


f (x) = x" — + .... 
= ax" — (n — 1) + =0 


By Sylvester’s theorem, the eliminant will be a relation of the form 
F (ki, ka, 


kn) = 0 (5) 


in which k,, ky, ....ky_, each occur to the degree n but k,, (or ¢) occurs 
to the degree n — 1, corresponding to the n — | stationary values. These 
(n— 1) stationary values will in general be different, but by a proper choice 


9) of ky, kg, ...-Akn, the eliminant (5) can be reduced to the form 
(¢ — A)" = 9, r+s=n—1 (6) 
Conversely, if (5) reduces to the form (6), the k’s are determined. The two 
™ stationary values given by (6) correspond to the ‘two extreme solutions’ 
referred to above. This offers a natural approach to our original problem 
but it seems difficult to carry out this programme in general. 
en 3. OTHER RELATED PROBLEMS 
The conditions for an extremum in the problem considered above can 
3) also be written as follows:— 
on 
(4) d . -Xn-) =0 
d(x;....Xn) =0 
m= the last of which holds because of the required stationary character of the 
function, while the remaining equations are consequences of the constancy 
of the concerned functions. It will be observed that exactly the same equa- 
tions are obtained if, instead of making ¢ = ky, stationary, any one of the other 
K's is made an extremum, all the remaining k’s being fixed. Thus the condi- 
“a tion (3) for an extremum holds for all the n different problems :— 


(i) any one k, = extremum, r = 1, 2,... orn 


(ii) all other ks, (n — 1) in number, (s+ r) fixed. 
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As in the original problem, we consider for each of the new problems only 
the case corresponding to the maximum number of double roots. By way 
of illustration we give the solutions for n = 3 and n= 4. 


The Case n = 3 


Problem a, 


k, =Stationary 
ky, ks fixed 


k,, = Stationary 
kg, ky fixed 


k, = Stationary 
ky, kz fixed 


Problem | a 


k, = Stationary m? (m—1) m? (m + 1) 


k, = Stationary (m?—m-+ 1)— (m* — m+1)+ 
ks, kg and k, fixed) (m+1)./m?+1 (m+ 1) 


k, = Stationary 1 m 1) 
ky, k, and k, fixed where a, c are the roots of nt he —2 yo 1) (m?—-1)t 
+ (3m — 1)? = 


k, = Stationary ead | 
ky, ky and k, fixed] m— | 


In our original problem (ky-stationary) we took a) = 0. It is readily 
verified that the Polynomial Identity formulated for the case a)#0 
transforms itself into the very identity for the case aj = 0 (p, p. 217, 220) 


by the mere substitution: 


dp — = a,’ 
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Taking y = 1, we see that, when a)+ 0, the differences (ay — a), (ay — a). -, 
(an — 4) are the roots of the Tschebyscheff polynomial found in p. We 
might here draw attention to one important distinction between the problem 
of ky*stationary and the others corresponding to k;(r#n) stationary. In 
the former case, the ratios (a, — ap): (a, — a): ....: (an — ag) are inde- 
pendent of any arbitrary parameter ; this is not the case for the other problems. 


4. REMARKS ON THE SYMMETRY PROPERTY 


For the solution of the original problem we established in p the following 
symmetry relations for any n: 


+ On = ay + = Gg + 


(8) 


It should be possible to prove these relations from the basic governing rela- 
tions, without obtaining the explicit form of the solution. 


Consider the case n = 2m +1. We have then the polynomial identity 
(p, p. 217), generalised to the present case ay+ 0, 


=(t — a9) (t — a9)? (t — a4)”. ..(t — oom)? 
= (t — a4)? (t — ag)?... .(¢ — Goma)? (t — + 


(€ — ao) (a, — a9)? (ag — — ay)*. 


Since this is an identity, a, a9,-a4,....a@gm as zeros Of (ft) are uniquely deter- 
mined as functions of a,, a3, ....¢2m—, §€. By writing € — t+ a, for ¢, the 
identity becomes, after a slight rearrangement 


(t — ag) [t — + — [t — + a — a5) ]* 
[t — (€ + a5 — Gem.) ]? 
= [t — (€ + a — [t — (€ + a9 — 
[t — (€ + a, — agm)]}*(t — 
+ (€ — a) (a3 — 


Comparing this with the previous form we can see that the identity is satisfied 
if we write 
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E + ay — = 


Of course, = agmi; = 


On account of the above-mentioned uniqueness, it follows that the 
sought relations between (ay, as, a4... and (a,, a3... are just 
these, which are of course the desired symmetry relations. However a similar 
proof does not work for the case n = 2m. 
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THE wave equation for a particle of arbitrary spin can be written in the form 
(Bhabha, 1949) 
+ X) =0, (1) 
where 
Pe=-—i x°(= x, x,? 


are the co-ordinates of a point in space-time; a are vx v matrices, and 
is a v-component wave function. X is a constant related to the masses of the 
particle described by (1). @° satisfies, in the general case, the minimal equation 


(a°)™ {(a)? — a,?} {(a°)? — a2}... — = 0. (2) 


Then also 


(P 


P™ (P2 — a,2p*) — a,%p%)... 


(3) 


ay? Pp 0, 
where 
P = Pk ak 
and 
= Px 
Equation (1) will describe a particle which has states of masses X/a,, 
Xi, 
The commutation relations for the particle field have been found by 
Udgaonkar (1953). They are 


(x), D] = —iS(x— x’), (4) 


where 


x Ar (x); (5) 
231 


S(x) = 


~ 
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here 4,(x) is the well-known invariant function of Jordan and Pauli for 
mass X/a,. [A, B] stands for AB + BA or AB—BA according as the 
charge density or the energy density is positive definite. 


We first prove that S(x) defined by (5) satisfies 
(P + X) S(x) = 


Now 


27.2 2 
(P+X)S(x) = > "A, P™ (P:— ) — x) 


where we have put 
Ar = 2 2 
(7) 
Since 4, (x) satisfies the equation 
(p? (x) =0, (8) 


we can replace X?/a,” by p? in equation (5); we get thereby 
(P + X) S (x) = 2A, P™ (P® — a,?p?) (P?— a,*p*)- -(P?— an? p?) A, (x) 


the last step follows from equation (3). 


Let 


A, (x)= — €(x) 4; (x), (9) 
where 
(x) = + 1 for x» > 1 


= — 1 for x5 < 1. 
Then [cf. Pauli, 1950] 


A, (x) = — 8(x), 


5 (x) = 8 8 (x) 8 (x*) 8 (x). 


L 
b 
r=1 
|_| 
2 
(10 
where 
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Let S (x) denote the function obtained from S (x) by replacing the 4, (x) 


by 4, (x), viz., 
x2 
S (x) = ZA, Pm (P*— (P — 
x(P p2 — Ay (x x) 
(11) 
We shall now prove that 
(P + X)S(x) = —(- 8(X). 
Operating on equation (11) from the left by (P + X) we get 
(P+ x)3(x) = ZA, Pm (p2— 4, (x) 
a) —SA,pm {pen (- 
—2) 2 1 pe 
x2 n 4 
-+(—1) \ A, (x). (12) 3 
Now from equation (3) : 
) pm (pan 4 pan-1)(_ g,2 p? — q,? p? — +++ — an? p?) 
‘ + (— 1)" p®” a,2a,? -an2} = 0. 
Substituting for P?"+™ from this equation in equation (12) we obtain 
= x2 x2 
(P+ x) S(x) =2A,;P™ [ + ag? =) 
x? 2(n—2) 27.2 xe 
+a,? + P {a, as on + 


ant {(ca) 4 


r 

| 
x. 
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On the right-hand side of this equation (p? — X?/a,) is a common factor. In 
view of the inhomogeneous equation (10) satisfied by 4, (x), we obtain 


(P + x) = — FA, P™| + a? + + ay?) 
x2 


2 x2 x2 \2 


Now the numbers A, defined by (7) satisfy the following identities [Bhabha, 
1950; cf. also the Appendix]: 


n 1 1 
2 Ar Ar = 9; G22) A, = (14 a) 


r=0 


(14 6) 


Hence on the right-hand side of equation (13) all terms, excepting the last 
one which contains (x?/a,”)"-', vanish. The contribution due to this term is 


m 
— (— 8 (x) 
so that | 


(P+ = —(— (15) 
We now define 


G(x) = [s (x) + x pm-2 + x2 


+(— (18 


+(— nm 3(x) |= — 


Hence G(x) is the Green’s function for equation (1). 


+ ( ) 1 2 n ar” ay” 
Then 
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The positive and negative frequency parts St (x) and S- (x) respectively, 
of S (x) are obtained by replacing the 4, in S(x) by 4,+ and 4,-, which are 


defined by 
sk 


i xt 


The function S“’ = i(S+ — S-) is also obtained by replacing the 4, in S by 
4,". All these functions satisfy the homogeneous wave equation. The 
Feynman function S, is defined by (cf. Pauli, 1950) 
Sr (x) = S® (x) — 21G (x) (17) 
and satisfies the equation : 


(P + X) Sp (x) = 21 8 (x). (18) 
Putting 
App = — 23 4, (19) 
we obtain 


r=1 


(ps — Aer (x) 
— — x2 PMs — (— 89), (20) 


To get the momentum space representation of S; (x) we make use of the 
following Fourier representations 


(21) 


8 (x) = dtp (22) 
px = p°x® — px! — px? — px3, 
p® = (p*)* — (p’)? — (p*)? — (p*)*. 


In the first equation above the integration over p', p*, p® is along the real 
axes, and over p® is along the real axis with two small detours one below 


= 
| 
Cc : 
2i 
Fr (27)* X2 
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the point — {(p’)® + (p?)? + (p*)? + X?/a,?}# and the other above the point 
+ {(p’)® + (p®)? + (p®)? + X2/a,2}t in the p® plane. We obtain 


Sr (x) = Ge d‘p (p) (23) 
where 


Sy (p) = -|z A,pm (P (pt — x) 


x x2 (Pe — x2 
+ 


(24) 
Making use of the identities (14) it can easily be proved that 


(P + x) Sy (P) = 1, 
so that 


Se(x) = f + (25) 


It should be noted that this is only a formal way of writing S; (x), since 
(P + X) is a singular operator for certain values of p;,, because (P + X) 4 = 0 
for non-vanishing ¥, and (P + X)-' then does not exist. At the poles of 


(P + xX) the path of integration in (25) has to be chosen as stated 
above. 


The function S; (x) given above is twice the propagation function K,. 
as defined in the theory of Feynman (1949). Hence the graphical techniques 
of Feynman and Dyson can be applied in the general case of particles of 


arbitrary spins if one uses the function 4 S, (x) for the internal particle 
lines. 


In conclusion, I should like to express my sincere thanks to Professor 
H. J. Bhabha for his interest in the work. 


: 
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APPENDIX 


The following proof of identities (14 a, b) is due to Dr. K. G. Rama- 
nathan : 


Let 
I(x) = (% — ay”) (x — (X — ay’) (A 1) 
and let ¢(x) be a polynomial of degree <n. Then 


(A 2) 


(A 3) 


(x) = x0? 


(ap? — ay”) (Ap? — (Gp? — Opa”) (Gp? — Gn?) 


(x — ay?) (x — ag?) — — Gp (X — 
(A 4) 


Equating the coefficient of x”-! on the right-hand side of this equation to zero 
we get 
(ar? — ay?) (ap? — ag?) — ya?) (Gp? — 
(A 5) 


This is the first of the identities (14 a) of the text. Similarly, taking 
= x”-3, x, we obtain other identities of set (14 a). 


To prove (145) take ¢(x)=1. Then 


1 
(a7 — Ars”) (a? — — an?) 


(x — ay?) (x — ag?) — apy?) — 


(x — (A 6) 


(x) Ay Ag A 

wher 

Ar=  (ar*) 
Let 

then 

n 2(n—2) 

1 = 

a 
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The constant term on the right-hand side must equal 1. Hence 


(A7) 
which is the identity (14 5). 
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INTRODUCTION 


IN a previous report (Lonappan, 1955) was described an investigation of 
the thermal expansion of potassium chlorate by an X-ray method. For 
this monoclinic crystal two of the principal expansion coefficients were small, 
of the order of 30 to 40x 10-8, while the third one was much larger, about 
130x10-*. The direction of the latter was almost perpendicular to the 
plane of the O; groups of the chlorate ions, all of which are parallel in the 
¢ structure of the crystal. It was thought worthwhile to extend the study to 
other crystals having similar O, groups, but having a different symmetry. 
Studies have been made with potassium nitrate (KNO,;) which is ortho- 
thombic and orthoboric acid (H;BO,) which is triclinic. The results with 
potassium nitrate are reported here; those with boric acid will be discussed 
in a separate paper. 


Potassium nitrate crystallizes in the class mmm and belongs to the space 
group — Pmcn. The unit cell dimensions are a = 5-43 A; b=9-17A; 
c=6-45A and there are four molecules per unit cell (Edwards, 1931). 
Since the crystal is orthorhombic the principal axes of the ellipsoid of expan- 
sion coincide with the crystallographic axes. The present study of thermal 
expansion has been made with a good degree of accuracy for a range of 
temperature 30 to 100°C. The study could not be extended for higher 
ranges of temperature, as in the case of potassium chlorate, since potas- 
sium nitrate undergoes a transformation at about 128°C., the crystal 
changing over to the rhombohedral form (Groth, 1910; Kracek, Barth and 
Kosanda, 1932) and the crystal, though retaining its external shape, gets 
broken into fine crystallites (Leonhardt and Borchert, 1936) as exhibited 
by the powder lines appearing in the rotation photograph above this tempera- 
ture. The complete ellipsoid of expansion could be determined by finding 
the variations in the spacings of the various reflections in [100] and [001] 
zones using appropriate crystals. 
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2. EXPERIMENTAL DETAILS 


Rotation photographs of the zero layer were taken for [100] and [001] 
zones using the 19 cm. camera (Unicam) and the spacings of the high angle 
reflections at the two temperatures were accurately determined using the 
85° fiducial mark in the camera, as described earlier (Lonappan, 1955), 
The overlap of reflections was avoided by the use of the ‘tilted crystal 
method’. The indices of the reflections determined by the reciprocal lattice 
diagram were confirmed by taking Weissenberg photographs of the crystal. 
The thermal expansion perpendicular to a number of crystallographic planes 
in a particular zone were thus determined. From these the principal expan- 
sion coefficients were calculated. Only reflections having a 9-value > 62° 
were used for the calculation of thermal expansion. Five such high-angle 
reflections were obtained in each zone. 


From a knowledge of the indices of the planes, the values of the inclina- 
tion ¢ of their normals to a chosen direction (b-axis) in the zone can be cal- 
culated. Then the relation between ag, the thermal expansion along any 
direction and the inclination ¢ of the direction concerned with b-axis takes 
the following form 


For the [100] zone, ag = ap + sind (ag — ap) (1) 
For the [001] zone, ag’ = ay + sin? d’ (ag — ap) (2) 


There are five equations of type (1) and five equations of type (2) from which 
@q, ap and ae are to be determined. Considering a as a function of ¢ and ¢’ 
three normal equations were formed using the ten equations, giving suitable 
weights to the a values according to their proximity to 90°. The probable 
errors were calculated using standard methods (Whittaker and Robinson, 
1944). 


As mentioned earlier, potassium nitrate undergoes a transformation 
at about 128°C. On heating single crystals of potassium nitrate above this 
temperature in the camera, the crystal became polycrystalline even when 
heated very slowly through the transition temperature and consequently the 
present method could not be used to study the thermal expansion for higher 
temperatures. The powder lines which appeared in the rotation photo- 
graphs about the [100] and [001] axes taken at 200° C. were measured and 
could be indexed on the basis of the rhombohedral unit cell of KNO, (1), 
a=4-50A; a= 73°30’ (Kracek, Barth and Kosanda, 1932). In both 
cases the crystals were found to give only powder lines when slowly cooled 
to room temperature and these corresponded to those given by the normal 
orthorhombic structure KNO,(II). The third modification of potassium 
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nitrate KNO, (III) reported to occur when KNO, is quenched from a high 
temperature was not observed in these experiments. 


3. RESULTS OF THE INVESTIGATION 


The measured values of the expansion coefficient perpendicular to various 
planes in the two zones studied are given in Table I. The principal expan- 


TABLE 


Measured thermal expansion coefficients of reflections in the 
[001] and [100] zones 


Indices Inclination Thermal expansion 
of reflection with b*-axis sin? d ap X 108 


660 59° 22’ —-0-7404 23-05 
1110 8° 44’ 0-0230 21-37 
570 50° 20’ 0-5925 23°19 
640 68° 27’ 0-8651 23-34 
2100 18° 40’ 0-1024 21-70 


038 75° 14’ 0-9350 171-73 
086 46° 50’ 0-5319 107-71 
028 80° 2’ 0-9700 177-85 
0 10 4 29° 38’ 0-2445 61-33 
095 38° 18’ 0-3841 84-52 


sion coefficients along the three crystallographic axes, calculated from these 
data, are 


ay, (b-axis) = 22:0 + 0°3; 
(a-axis) = 23:5 +0°5; 
agg (c-axis) = 182-6 + 0:2; x 10-6 
Fig. 1 brings out the agreement between the observed values and those cal- 


culated from the constants of the ellipsoid of expansion. In Fig. 1, a is 
plotted against sin? ¢ or sin? ¢’ as the case may be. The observed values 
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1604 


Fic. 1. Variation of expansion coefficient (30 to 100° C.) with sin? @ in the two zones [100] 
and [001]. The circles represent the experimental values and the straight lines have been drawn 
using the calculated values of a,, a, and a,. 


should lie on two straight lines and they are found to be very close to the 
theoretical values. 


It is interesting that ag and ap are found to be almost equal whereas a, 
is very much higher. This is due to the fact that the planes of all the nitrate 
ions occurring in the crystal are parallel to one another. According to the 
structure determined by Edwards, the O, groups all lie in the c-plane. In the 
present study the maximum value of thermal expansion coefficient is obtained 
along the c-axis, which is perpendicular to the plane of the O, group, as it 
should be. The similarity between the morphology and structure of potas- 
sium nitrate and aragonite has already been pointed out by Miller (1840) 
and Edwards (1931). However the actual values of the thermal expansion 
coefficients for aragonite are 10, 16 and 33 x 10-® while for KNO, they are 
22, 23-5 and 183 10-*. It will be noticed that the value of the largest co- 
efficient is very different in the two cases. This is not surprising for a similar 
large difference occurs even between calcite and NaNO, (Table II) which 
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are closely related in structure. It appears that the coefficient normal to 
the plane of the ions is much larger with nitrates than with carbonates. 


TABLE II 
Thermal expansion of some crystals of the type ABO, 


Thermal expansion x 10° 
Crystal 

digg * 
Aragonite ia ea 10 16 33 
Calcite .. we a — 6 — 6 25 
NaNO, .. A 11 120 
KCIO, .. 30 39 134 
KNO, .. 22 183 


*as, always normal io plane of O, group. 


Coefficient of cubical expansion calculated from the present data namely 


Geuy = (a4; + aa + agg) = 228X10-° agrees with the directly measured 
value of 21110-® reported by Bellati and Finazzi (1910). 


The value of cubical expansion coefficient may be used to verify the 
Griineisen relation for potassium nitrate. The density of the crystal at room 
temperature is 2-098, from which the molar volume V, is calculated to be 
48-2. The molar specific heat (C) as given in the Landolt-Bornstein Tabellen 
is 24-15. Unfortunately the value of the compressibility K, required in the 
relation 

¥ = Vo%cur/KoCv 


(where y is the Griineisen constant) is not available in the literature. How- 
ever y may be evaluated from the formula 


1 dv 
_ dog») _ _ iv iX, 
VoT 


Here X, is the proportional rate of change of frequency with temperature, 
as defined by Ramachandran (1947 a). The value of X, for potassium nitrate 
is not available, but the corresponding value of sodium nitrate may be cal- 
culated from the Raman Effect data for this crystal (Nedungadi, 1939). 
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The tilting and transverse oscillations of NO, group in sodium nitrate have 
X-values of 800 and 600% 10-® respectively (Ramachandran, 19475), the 
mean of which, namely 7x10-* may be taken for our calculations. This 
yields the value 3-07 for the Griineisen constant. 


Putting this value of y in the equation mentioned earlier, we get 
K, = 35-28x10"* 


This value for compressibility is of the same order as for other nitrates, 
namely 31-2 for barium nitrate, 34-2 for strontium nitrate, 28-2 x 10-'S for 
lead nitrate, which indicates the approximate validity of Griineisen’s rela- 
tion for the crystal. 


SUMMARY 


The principal coefficients of expansion have been determined for the 
orthorhombic crystal potassium nitrate for the range of temperature 30 to 
100° C. by an X-ray method. The values are a,; (b) = 22:0; ag. (a) = 23-5; 
agg = 182-6 10-*. The direction of the greatest expansion is per- 
pendicular to the plane of the O, group and the expansion is nearly constant 
for all directions in the plane of the O,; group. The approximate validity of 
the Griineisen relation is tested by calculating the compressibility from a 
knowledge of the variation of the Raman frequency with temperature. 


The author wishes to record his deep indebtedness to Professor G. N. 
Ramachandran for the invaluable suggestions he gave and the keen interest 
he took in the course of this study. 
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STUDY OF THE ANISOTROPY OF COSMIC RAYS 
WITH NARROW ANGLE TELESCOPES 
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(Physical Research Laboratory, Ahmedabad-9, India) 
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WHEN appropriate corrections to remove terrestrial influences are applied 
to the daily variation of meson intensity, it shows up the anisotropy of the 
primary cosmic radiations. as viewed by an apparatus fixed to the spinning 
earth. Early data relate almost exclusively to measurements made with 
ionisation chambers for study of the various types of time variations of cos- 
mic ray intensity. In these the instruments measure radiation from all direc- 
tions and are therefore not particularly appropriate for study of the daily 
variation of intensity. Several later studies have been made with counter 
telescopes having directional characteristics. However in an effort to secure 
a high counting rate and good statistical reliability of the results, the tele- 
scopes generally subtend large angles and hence most of them ignore the 
requirements for a specialised study of the daily variation. From available 
determinations it has been in consequence erroneously assumed by several 
workers that the meson intensity at sea level does not exhibit a daily variation 
with an amplitude in excess of 0-2 to 0-4%. The notable exceptions to this 
way of thinking have been the Japanese group of workers and the group at 
Physical Research Laboratory, Ahmedabad. 


Sekido and his co-workers have drawn attention in a series of papers'?* 
to the remarkable differences observed in the nature of the daily variation 
and its changes during periods of magnetic disturbance when measurements 
are made with vertical telescopes having semi-angles of 12°, 40° and 85°. 
Having concluded that the daily variation is better observed with narrow 
angle telescopes they* have conducted an experiment with a telescope having 
a semi-angle of 5° to determine sources of cosmic rays in the galaxy. 


About three years ago experiments were commenced at Ahmedabad to 
determine the result of pushing the technique of narrow angle telescopes to 
what was considered a practical lower limit. Observations have been made 
by one of the authors (P. D. B.) with a triple coincidence telescope having 
a semi-angle of 1-8° in the E-W. plane and a semi-angle of 6-7° in the N-S. 
plane. The intensity relates to the component of cosmic radiation which 
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can penetrate a minimum thickness of about 27 cm. of iron at sea level. The 
experiment was conducted from November 1952 to August 1953 and the 
mean bi-hourly rate was about 35 counts. The percentage mean daily varia- 
tion is shown in Fig. 1. In view of the large probable error of each bi- 
hourly deviation, moving averages over three consecutive bi-hourly values 


1.8 


00 04 08 12 16 20 00 
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Fic. 1. Daily variation measured with 18T. 


have been taken. Harmonic analysis of the daily variation given by the 
bi-hourly percentage deviations from mean, before taking moving averages, 
shows that the diurnal amplitude M” is 2-2 + 0-5% and occurs at noon. 
The semidiurnal component with an amplitude of 0-4 + 0-5% is not signifi- 
cant but has a negative correlation of — 0-81 with the semi-diurnal compo- 
nent of the daily variation of barometric pressure. 


The large amplitude of the daily variation found in this experiment with 
a narrow angle telescope has prompted us to make a more elaborate study 
involving simultaneous measurements with telescopes of different angles. 
It was hoped thereby to determine with some precision the profile of the 
anisotropy of primary cosmic radiation and also to establish an optimum 
experimental technique for study of the daily variation of meson intensity. 


In Fig. 2 we show the arrangement of counters in an experiment set up 
by one of the authors (N. W.N.). The apparatus furnishes triple coincidences 
from three independent telescopes 2°T having semi-angles of 2-5°, from two 
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Fic. 2. Experimenta! Arrangement. 
5T —(A’B’) (A’B’) (Ar (B’C’) (B’C’) 
bT —(A’B’C’”’) (A’B’C’) (A 


telescopes °T with semi-angles of 5° and from one telescope “T with a semi- 
angle of 15° in the E-W. plane. All telescopes have semi-angles of 19° in 
the N.-S. plane. In addition, the total rate °T from a battery of counters 
has been taken to give the omni-directional intensity. 12cm. of lead has 
been used as shielding in every case. To increase the counting rate of tele- 
scopes without change of semi-angle, duplicate sets of counters are placed 
horizontally displaced from the first set, but connected with them in parallel. 
This introduces no error except for a negligible contribution from penetrating 
air showers. 


In Table I are given details of the telescopes, their characteristics and 
period of operation. Since it is known that the daily variation of meson 


TABLE I 
Details of Telescopes 


Telescope 247 
Daily Variation .. 25M 5M 16M 90M 


Semi-angle of Telescopes... 2°5° 15° 90° 
No. of Telescopes me 6 4 1 1 
Period for which data availa- 1954 1954 1954 1954 
able Jan.-May. Jan.-May Sep-Dec. Jan.-May 
Sep.-Dec. Sep-.Dec. Sep.-Dec. 
1955 1955 1955 1955 


Jan.-March Jan.-March Jan.-March Jan.-March 


Mean bi-hourly counting rate 
per Telescope .. 199 251 x4 448 x 16 499 x 1024 
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intensity can alter with the passage of time, it is obviously necessary to com- 
pare the performance of the different telescopes over a period when all 
instruments have simultaneously been working. For such a period we show 
in Fig. 3 the mean percentage daily variations *°M, °M, °M and %M, 
superposing the data of as many similar telescopes as have been simulta- 
neously in operation. Figures 3a and 35 relate respectively to the daily 
variations before and after application of a barometric pressure coefficient 
Bi —2-2% per cm. of Hg. The amplitudes and the times of maxima 
of the diurnal and the semi-diurnal harmonic components of the daily varia- 
tions measured with different arrangements are given in Table II. The diurnal 
components are also shown on harmonic dials in Fig. 4. The values before 
and after correction for the daily variation of barometric pressure are shown 
alongside of each other, the former being indicated by asterisks. 


The following features emerge from a comparison of the daily variations 
of meson intensity observed with the different types of instruments. 
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Fic. 3. Daily variations measured with *°T, 5T, and T. Figures (a) and (5) relate 
respectively to values before and after applying correction for daily variation of pressur?. 
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Fic. 4. Harmonic dials showing diurnal components of 25M, 5M, 15M and °M before 
(L.H.S.) and after (R.H.S.) applying correction for daily variation of pressure. 


1. *5M and °M are almost identical. They differ markedly from “M 
and again from 9M. 


2. The percentage amplitude of the diurnal component. *°M? is slightly 
greater than °M”, but not significantly so. However the amplitude which 
is about 1-2% for the narrow angle telescopes, decreases rapidly for wider 
semi-angles beyond 5°, and is only about 0-25% for the omnidirectional 
intensity. 

3. The time of maximum M¢?” of the diurnal component corresponds 
to about 0900 hours for **°M and *M but becomes progressively earlier with 
wider semi-angles. °*M¢? corresponds to about 0500 hours. 


4. The semi-diurnal component of the daily variation does not alter 
as much as the diurnal component when measurements are made with the 
different instruments. There is an indication that the amplitude M§ decreases 
with increasing angle of telescope, but not to the same extent as M?. The 
time of maximum is fairly constant. 


5. There is high negative correlation between the semi-diurnal varia- 
tions of meson intensity and of barometric pressure, the omnidirectional 
intensity having the highest correlation. The application of a uniform 
correction for the daily variation of barometric pressure using a coefficient 
B = —2-2% per cm. of Hg, results in a virtual elimination of the semi- 
diurnal components in 7°M, °M and ™M. However %M is still left with a 
significant semi-diurnal component. 


The choice of 8 = — 2-2% per cm. of Hg is dependent on considera- 
tions described elsewhere,® but it is not clear whether the application of the 
same coefficient to °°T as well as to the narrower telescopes is fully justified. 
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It will be observed that in *°T and °T, the barometric correction with a co- 
efficient of this magnitude makes relatively little difference to the daily variation 
of meson intensity. Hence, there is likely to be no great ambiguity in inter- 
pretation on account of the operation of uncertain meteorological influences. 


Computation shows that meson intensity incident on °°T from direc- 
tions inclined less than 5° from the vertical, makes a contribution of less than 
1-0% to the total omnidirectional counting rate. Thus %M with a diurnal 
amplitude of 0-24% represents the integrated daily variation of intensity 
incident from all directions, and is little affected by the special features of 
the daily variation characteristic of intensity in directions close to the vertical. 
On the other hand 25M and °M indicate that the intensity which exhibits a 
large daily variation of diurnal amplitude | -2% is confined almost exclusively 
to a narrow vertical cone of incidence with a semi-angle of about 5°. 


If we assume that intrinsically there are two types of daily variations, 
the V-type characteristic of the vertical intensity and the O-type characteristic 
of the omnidirectional intensity, then. we have V ~ 5M — °°M and O°°M. 
These are shown in Fig. 5. 'M is compatible with a superposition of the 
two types with appropriate weightage factors. This is clear from Figs. 3 


V 


+ 1-0] 


6 12 18 
HOURS 1.S.T. 
Fic. 5. The V and the O types of daily variations of meson intensity. 
and 4 where we show the experimental and computed values of *M and its 
diurnal component respectively. The weightage factors have been evaluated 
from considerations of the geometry and dimensions of °T, T and °T, 
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It is remarkable that there is sharp attenuation of intensity showing 
V-type of daily variation for incidence at angles exceeding 5° with the vertical. 
Simpson et al. have reported for local production of neutrons a daily varia- 
tion comparable in amplitude to V. Since the nucleonic component of cos- 
mic radiation, which is responsible for the local production of neutrons, gets 
rapidly attenuated with increasing zenith angle, the neutron monitor essen- 
tially measured vertical intensity. Thus a large amplitude of the daily varia- 
tion appears to be a characteristic of a radiation, which can only penetrate 
to low levels of the atmosphere at small zenith angles. The main difference 
in amplitudes of the daily variations* observed with Simpson’s neutron 
monitor and the Sitkus ionisation chamber could possibly arise from con- 
siderations of directional sensitivity of the measuring instruments rather 
than from difference in mean energies of primaries to which they respond. 


It is important to understand the O variation and its origin. This is 
dependent on study with narrow angle telescopes pointing towards directions 
making different angles with zenith in the E.-W. plane. If, as appears likely, 
the O variation does not arise due to an integration of daily variations vastly 
different in character for different zenith angles towards east and west, we 
might be here dealing with a general modulation of the cosmic ray intensity. 
This would presumably occur near the top of the atmosphere but be under 
solar control. 


Narrow angle telescopes open up a wide new field for experimental 
investigation of the anisotropy of primary cosmic rays. With low counting 
rates, data have to be averaged over an extended period to get significant 
results. Since the nature of the anisotropy changes, often radically, with 
passage of time, the averaging of the daily variations for several days is rather 
unsatisfactory. We can however, get around this difficulty by having simul- 
taneous measurements with an array of several independent narrow angle 
telescopes. Work is now in progress along these lines and a report will be 
published separately. 

We have benefited by discussions with Dr. R. P. Kane. We gratefully 
acknowledge support of the Atomic Energy Commission of India to our 
project. 
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THE OCCURRENCE OF HIGH MULTIPLE 
REFLECTIONS FROM THE F2 REGION OF THE 
IONOSPHERE BASED ON A STUDY OF THE 
AHMEDABAD RECORDS 


By R. G. Rastoai, M.Sc. 
(Physical Research Laboratory, Ahmedabad) 
(Communicated by Dr. K. R. Ramanathan, F.a.sc.) 


INTRODUCTION 


SHorTLY after the installation of the automatic ionospheric recorder at 
Ahmedabad in January 1953, it was observed that on certain nights groups 
of pulses other than the regular ones moved swiftly across the time base 
on the monitor oscilloscope. On the P’-f records, they produced traces 
which started from below the ground pulse and crossed the various P’-/ traces. 
Examples of such records are shown in Fig. 1. After the first three months 
of regular P’-f recording, the phenomenon was not observed for some months, 
but reappeared clearly in December 1953 and January 1954. A study of 
the phenomenon was then undertaken. 


GENERAL FEATURES OF THE PHENOMENON 


The following features were noticed about these abnormal traces: 
(1) They occurred only during night. (2) At any frequency the separation 
in height between two adjacent abnormal traces appeared to be the same 
as the virtual height of the first reflection trace at that frequency. (3) They 
ended at the ordinary critical frequency f,F, of the F, layer, and on some 
occasions, there was another group of pulses ending at the extraordinary 
critical frequency fF, of the F, layer. (4) The ratio of the slope of an 
abnormal trace to that of the first reflection trace was found to be the same 
at any frequency and increased in steps of one from one trace to the next. 
It was therefore concluded that these abnormal traces represented high order 
reflections from the ionosphere. 


The equipment at Ahmedabad is a British N.P.L. type automatic 
ionospheric recorder and the pulse repetition frequency is 50 per second. 
If the delay time of a high order multiple reflection is more than 1/50 second, 
it will be recorded in the second sweep of the time base, and the equivalent 
height of the multiply reflected trace will exceed 3000 km. (= 1/50 second). 
The order of reflection can be determined by dividing the slope or the equiva- 
lent height of the trace at any frequency by the corresponding slope or height 


_Of the first reflection trace at that frequency. 
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During regular recording, the duration of the time base is usually 
adjusted to record echoes of equivalent height upto 1000 km. only. How- 
ever, the time base duration can be extended to record echoes of equivalent 
height upto about 2000km. With such an extended time base, multiples 
upto the 17th order have been recorded. It was considered that it would be 
worthwhile to increase the height range still further and arrangements were 
made on a few occasions to record the echoes in two steps, first showing the 
reflections upto 1500 km., and then delaying the start of the time base so as 
to record heights between 1500 and 3000km. Such a double record is 


shown in Fig. 2 where all multiple echoes from the first to the ninth can be 
seen. 


An examination of the literature showed that high multiple reflections 
have been noticed by some previous workers on P’-t records (Pierce and 
Mimno,! Baird’). They attributed these to the curved shape of the iono- 
sphere or humped ionisation contours in the ionosphere. These humps 
produce a mirror-like converging action on the reflected radio beam. Gipps, 
Gipps and Venton® also explained some abnormal traces in P’-f records in 
terms of humps and concavities in the ionosphere. Munro‘ described certain 
complexities observed on ionospheric records as the effects of the curvature 
of the ionospheric surface. Uyeda and Nakata® have recently described 
similar abnormalities observed on h’-t and fct records and explained them 
on an assumed inclination of the ionosphere and its movement. 


ANALYSIS OF AHMEDABAD RECORDS SHOWING HIGH ORDER 
MULTIPLE REFLECTIONS 


A detailed analysis was made of the records showing such abnormal 
traces and the present note contains a summary of the findings. 


The parameters of the F, layers were determined by the method of 
Appleton and Beynon.* The virtual heights of the layer at various fre- 
quencies were plotted against the corresponding values of Booker and Seaton’s 


function’ 
( = f, Be 


where f; is the critical frequency of the layer and f any frequency less than 
the critical frequency. It was found that a neat straight line could be drawn 
through the points showing that the distribution of electrons in the F, iayer 
during the night time was approximately parabolic. The values of the virtual 
heights at ¢ = — | and ¢ = O gave respectively the true height of the base 
hoF, and the height of the level of maximum ion density hpF,., the semi-thick- 
ness being given by Ap— My. 
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The analysis of a typical record (Fig. 1 5) is given in the following table. 
The value of f,F. was 3:0 Mc/s. h(n), h(n), h(ns) and h(n,) are the 
virtual heights of observed multiples and are found to be nearly equal to 
11, 12, 13 and 14 times 4,. ApF, was found to be 250 km., AgF, to be equal 
to 198 km. and the semi-thickness of the layer was 52 km. 


With a view to find out what ionospheric features were associated with 
high multiples, the parameters ho, hp and T obtained from the hourly iono- 
spheric records during the occurrence of the phenomenon were examined. 


30 


Total Number of Occurrences 


(b) 
= 
£ 20F 

260 
(a) 4 


> ae 20 22 00 02 04 06 
75° E.M.T. in hr. 


Fic. 3. (a) Mean minimum virtual height of the F,-layer (January 1953 to December 1954). 
(b) Average hour to hou: change in h’F, (January 1953 to December 1954). (c) Total number of 
occurrences of high multiples (January 1953 to December 1954). 
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In the records of the period January 1953 to December 1954 there were 256 
such occurrences. It was found that in general high multiple reflections were 
associated with a hump in the ionospheric surface over the place of observa- 
tion. High multiples occurred more or less continuously for many hours 
in certain periods of the night in some months. Near the critical frequency 
the multiples are generally weaker, due no doubt to increased absorption 
as well as to increased slope of the P’-f traces. 


The distribution of the occurrence of high multiples with the hour of 
the night is shown in Fig. 3 (c). In Fig. 3 (a) and 3 (b) are shown the mean 
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Fic. 4. Variation of hgF, and the occurrence of high multiples on a few nights. 
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virtual height (h’F,) and the average hour-to-hour change in virtual height 
(h'F,) during the period January 1953 to December 1954. 


It will be seen that the hour of maximum occurrence of high multiples 
are round about 21 hours, 01 hour and 04 hour. The times of maximum 
frequency of occurrence of high multiples nearly coincide with or precede 
the times of maximum rate of change of h’F,. The strength of the peak at 
01 hour in Fig. 3(c) is abnormal due to the large number of occurrences 
of high multiples at that hour in one month only, namely in November 1954, 

Figure 4 gives the analysis of the heights of the F, layer on a number 
of individual nights at times when high multiples were observed. The circles 
represent the times at which high multiples were observed. 


Observations at short period intervals taken on a few nights showed 
fair correlation between the highest order multiple and the curvature of the 
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Fic. 5. The variation of the total number of multiple reflections and the height of base 
of the F, layer on the night of 10 to 11 March 1955. 
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ionospheric surface. Figure 5 shows the results of observations taken on 
the night of 10-11 March 1955. 

It was found that sporadic E does not affect the high multiples except 
when it blankets the normal reflections. Sometimes high multiples of F, 
and sporadic E can be seen simultaneously recorded over the same frequency 
range (Fig. 15). 

1 F scatter usually reduces the number of multiple reflections but there 
are some exceptions to this also (Fig. 1 c). 


of Occurrence 


Percentage Frequency 


SEASONAL VARIATION OF MULTIPLE REFLECTIONS 


It appears from Fig. 6 that high multiples are least frequent in the 
summer months May to July and most frequent in November to March or 
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Fic. 6. Percentage of nights on which high multiple reflections were recorded in various 


months January 1953 to December 1954. 
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April. This accords with the fact that the daily variations of h’F, are least 
pronounced in the summer months and most pronounced in autumn and 
spring. During summer, E; and 1 F—scatter are rather more frequent and 
these no doubt contribute partly to the rarity of high multiples in that season. 


IRREGULARITIES IN THE INTENSITIES OF THE MULTIPLE REFLECTIONS 


A phenomenon related to the above is the anomalous and _ varying 
distribution of intensity of multiple echoes. This phenomenon can be seen 
clearly on many occasions on the monitor oscilloscope and generally occurs 
on multiples higher than the third. It is difficult to determine quantitatively 
intensity ratios of the various multiples from the P’-f records but approximate 
estimates of their intensity can be obtained in some cases. 

At times, depending upon the position of the hump with respect to the 
place of observation, a higher multiple may be formed by combination of 
waves leaving the transmitter at different angles and arriving at the receiver 
after being reflected from off-vertical points in the ionosphere. The intensity 
of some particular multiple may then be greater than that of the lower ones 
depending upon the curvature, height, and reflection coefficient of the F, layer, 
as well as on the radiation patterns of the transmitting and receiving antenna. 

CONCLUSION 

These observations show clearly that even at night, the F, layer can 
by no means be considered as a simple plane reflecting surface and that there 
are dynamic changes going on in it most of the time. Determination of the 
reflection coefficient of the layer and of the absorption of the atmosphere 
below it is rendered particularly difficult because of these changes. 
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FIG, 1(4). P’-f record at 02 hr. on 31-1-1953, 


FIG, I(c). P’-frecord at 00 hr, on 17-3-1953. 


P’-f records showing abnormal traces crossing the lower order reflection traces. 
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P’-f record at 22 hr. on 10-3-1955 recorded in two steps, one recording echoes upto 


and the second immediately after, from 1500 km. to 3000 km. 
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FIG. (a) P’-f record at 00 hr. on 18-3-1954. 
(2) 
FiG. (4) Photo at 2230 hr. on 28-3-1955, 
FiG. 8. Records showing second order reflection stronger than the first order reflection 
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(2) 
2125 hrs. 


(4) 
2126 hrs. 


(c) 
2127 hrs. 


(2) 
2128 hrs. 


FiG. 9. Multiple echoes at intervais of one minute on 24-2-1955. 
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